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Abstract. The detection of impulsed TeV emission from 
PSR B 1706-44 raised the question if this VHE emission re- 
sults from inverse Compton scattering of multi-TeV elec- 
trons on the CMBR and other soft photon fields. This 
hints at the existence of an unseen synchrotron nebula, 
which is larger than the observed V compact X-ray neb- 
ula. In a search for the synchrotron counterpart of the 
TeV nebula, we have taken V-band CCD images of the 
region around the PSR B 1706-44 pulsar position. By re- 
jecting all point source contributions down to the diffuse 
sky limit, we were able to search for such extended ple- 
rionic synchrotron emission on scale sizes limited by the 
angular resolution of the CANGAROO imaging TeV ob- 
servations. By combining the V-band flux limits with the 
observed TeV flux, we were able to obtain upper limits for 
the magnetic field strength as a function of the radius of 
the assumed TeV nebula. Assuming that the 2' radio ple- 
rion also defines the size of the TeV plerion, we constrain 
the steady state particle spectral index to values between 
4 and 5.5, with a field strength limit of about 3 fiG. This 
compares with an index of 3.6 ± 0.6 inside the V compact 
X-ray nebula. 

Key words: Magnetic fields - Methods: data analysis 
- pulsars: individual: PSR B 1706-44 - Radiation mecha- 
nisms: non- thermal - X-rays - gamma-rays 



1. Introduction 

PSR B 1706-44 is a young Vela-like pulsar (spin-down age 
17.5 kyr), with a period of 102 ms and a large spin-down 
power (E = 3.4 x 10 36 ergs/s). It was discovered during a 
20 cm radio pulsar survey of the southern Galactic plane 
(Johnston et al. 1992), and later detected in soft X-rays 
during the ROSAT mission (Becker et al. 1992, 1995) and 
identified as a pulsed GeV source by EGRET (Thompson 
et al. 1992). Very high energy (VHE) 7-ray observations 
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above 1 TeV from CANGAROO (Kifune et al. 1995; Ki- 
fune 1997) and above 0.3 TeV from Durham (Chadwick et 
al. 1997) detected and confirmed the existence of unpulsed 
radiation of statistical significance from this source. The 
CANGAROO detection is consistent with a point source, 
with angular extend not exceeding ~ 0.°12 - the pixel 
size of the imaging camera. A dispersion based distance 
measure of Taylor & Cordes (1993) places PSR B 1706-44 
at ~1.8 kpc. A possible SNR G 343.1-2.3 association with 
the pulsar was proposed by McAdam et al. (1993), but 
this was later found to be unlikely by Frail et al. (1994) 
and Nicastro et al. (1996). 

Chakrabarty & Kaspi (1998) gave a red-band upper 
limit of R > 18 to the pulsar. A 3<r upper limit to the pul- 
sar magnitude of V = 24.5 was given by Lundqvist et al. 
(1999), and this is consistent with the theoretical predic- 
tion of V = 24.12 (Urama & Okeke 1998, and references 
therein). Mignani et al. (1999), using the same data as 
Lundqvist et al., obtained an upper limit of V > 27.5. 

The photon flux above 1 TeV from this unpulsed source 
is only two times smaller compared to the flux of the Crab 
Nebula at the same energy. This is remarkable in the con- 
text of a synchrotron-inverse Compton (on external pho- 
ton fields) interpretation for the TeV 7-rays from PSR 
B1706-44, since the synchrotron nebula of this source is 
very weak compared to the Crab Nebula's synchrotron in- 
tensity. Whereas the Crab Nebula's intense synchrotron 
emission is the result of a large B (~ 10 -4 G) field, we 
must have a much weaker field for this plerion to avoid a 
bright synchrotron nebula. 

De Jager (1995) speculated that a Vela-like compact 
synchrotron nebula (< 1', Harnden et al. 1985) may be 
present, which would account for most of the unpulsed 
X-ray emission from this source, since the pulsar wind 
magnetic field (which scales as B oc 1 / R) would predict a 
compact synchrotron nebula, if the conditions are similar 
to Vela. The corresponding particle density required for 
the synchrotron nebula would however be too low to pro- 
duce a detectable inverse Compton compact nebula, since 
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the target photon density from the CMBR and Galactic 
disk would be too small. 

A compact nebula was indeed discovered (Finley et al. 
1998), which confirms the abovementioned interpretation. 
The scenario of Harding & de Jager (1997), with a de- 
tailed analysis by Aharonian et al. (1997) may apply in 
this case: electrons are streaming away from the compact 
X-ray nebula into a low-B extended plerion. If the dif- 
fusion coefficient for electrons in the extended plerion is 
small enough, it may be possible to trap enough electrons 
to account for the TeV emission by IC scattering. 

Frail et al. (1994) detected a radio synchrotron nebula 
associated with the pulsar, which has a radius of ~ 2' 
at 20 cm wavelength. This may represent the low energy 
counterpart of the TeV nebula, since the associated B in 
the radio nebula is smaller than the pulsar driven field in 
the compact X-ray nebula, and the size is just smaller than 
the angular resolution of the CANGAROO telescope. 

If the TeV 7-rays originate from the extended (radio?) 
nebula, we will expect a bright extended nebula in optical 
and/or X-rays, if the field strength is large enough. No ex- 
tended X-ray nebula was however seen, resulting in an up- 
per limit of four times the flux of the point source/compact 
X-ray nebula at 1 keV, for R = 2' (K. Brazier 1998, per- 
sonal communication). This is consistent with the expec- 
tation that the pulsar field strength (which drops as 1/R 
outside the light cylinder) should start to drop below the 
ambient field strength of a few fiG at a distance of a 
few arcminutes (de Jager & Harding 1998). The conse- 
quence of this is that both the synchrotron brightness and 
synchrotron characteristic frequencies should drop with 
increasing R. Scaling from the Durham detection above 
0.3 TeV, and the expected B ~ 3 fiG in the extended 
nebula (de Jager & Harding 1998), we may use equation 
(5) of Aharonian et al. (1997) to calculate the synchrotron 
frequency which corresponds to the Durham detected 7- 
rays near 0.3 TeV, assuming that the CMBR is mostly 
responsible for the inverse Compton scattering. This gives 
a frequency of 

v ~ 1.5 x 10 15 ( El \ ( Hz, 

which is within a factor of 3 from the V-band. The non- 
detection of the unpulsed component at 20 GeV (Thomp- 
son et al. 1996) hints at a marginal turnover in the 7-ray 
spectrum (de Jager & Harding 1998), which corresponds 
to a turnover frequency as low as 10 14 Hz, which is well 
below the V-band. 

At a distance of 1.8 kpc to the source, we expect sig- 
nificant interstellar absorption in the blue and UV bands. 
We therefore take wide field V-band images as a compro- 
mise between interstellar absorption and avoiding a spec- 
tral turnover at too low frequencies. By taking overlapping 
images, we cover a total field-of-view of about 10', which 
should include the expected TeV source. Even if the full 
extend of the TeV source is 10', we still expect a typically 



Table 1. Observing log for B 1706-44. The 1st frame is 
the central one. 



Frame RA (J2000) Dec MJD (start) K time(sec) 



1 


17 09 42 


-44 28 57 


51322.80399 


29.7 


720 


2 


17 09 52 


-44 31 20 


51322.81410 


28.7 


720 


3 


17 09 33 


-44 26 34 


51322.82833 


29.2 


720 


4 


17 09 52 


-44 26 34 


51322.83872 


28.3 


720 


5 


17 09 33 


-44 31 20 


51322.85028 


28.6 


720 



centrally brightened image, resulting in a radial gradient 
from the source. Failure to detect the optical counterpart 
will allow us to set upper limits on the magnetic field 
strength in the extended nebula, since any stronger field 
should have resulted in a brighter optical nebula. 

2. Observations and data analyses 

V band CCD observations of PSR B1706-44 were carried 
out on May 24, 1999 at the South African Astronomical 
Observatory (SAAO) 1.0 m telescope in Sutherland. Five 
frames, each ~ 5. 3 x 5. 3 (pixel size = 0.31 arcsec), were 
taken. The first one was centered on the radio pulsar po- 
sition (Johnston et al. 1992, 1995), whereas the other four 
frames were offset in such a way that they have an overlap 
of ~10% with each other around the pulsar position. The 
overall region observed around the pulsar is 9. 5x9. 5. The 
observation time for each frame was 12 minutes. Table [l] 
gives the position (RA, Dec), the starting time (MJD), 
K values and the duration time (seconds) log for the five 
frames. 

Observations were affected by the moonlight (about 
60% bright). The average counts for each frame were de- 
creasing with time, probably due to the moon's effect 
across the sky. The data were cleaned and flat-fielded us- 
ing the IRAF image-processing software. Fig. |l| shows a V 
band CCD image of the central frame with the pulsar and 
the standard star positions labeled + and 1 (Mignani et 
al. 1999, their Figs. 1 & 2), respectively. Only the inner 
part (25% of the total area) of Frame 1 is shown on Fig. [l| 

For each frame, we made a histogram of the number of 
pixels versus intensity per pixel using a logarithmic bin- 
ning procedure prescribed by 

K = N \og{N/N rnin ), 

where N Q (= 99/ \og(N rnax /N rnin )) is a normalizing con- 
stant, as shown in Fig. || for the central frame. This bin- 
ning procedure ensures that all intensities are confined to 
the interval K = 1 and 100. 

Whereas the diffuse emission is expected to be confined 
to a rather narrow low- if peak in the histogram, the pixels 
contaminated by point sources should show a broad tail. 
Fig. U indeed shows such behaviour. It can be shown that 
the weak plerionic emission is expected to contribute to 
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Fig. 2. A non-linear intensity distribution histogram of 
the central frame (example). K = 1 corresponds to 264 
counts (minimum), whereas K = 100 corresponds to 
Fig. 1. Star 1 is a standard star, + is a pulsar position 6.89 x 10 4 counts (maximum). The cut is made at K ~ 30, 



(Mignani et al. 1999). The frame is zoomed in. The origi- 
nal frame is about 1024 x 1024, CCD frame size. 



above which the high intensity star light contribution is 
rejected. 



the low K- values (K < 25), so that we can reject the star- 
like objects by cutting at K > 30 (for frame 1) as shown 
in Fig. [|. This rejects only 10% of the pixels, so that we 
preserve 90% of the plerionic component. The threshold 
value of K for each frame is given in Table [lj and proba- 
bly depends on the phase of the moon and the background 
bright stars (in each frame), which contribute to the back- 
ground for each frame. Note that the final results are not 
sensitive to the threshold K. 

By renormalising the count rate of each frame to the 
count rate of the central frame, we were able to construct 
a composite image of a ~ 10 xlO region around the pul- 
sar. Since the plerionic structure is expected to be circular 
around the pulsar (the radio nebula is at least circular), 
we constructed annuli around the pulsar and made a plot 
of the mean count rate per pixel versus radius of the an- 
nuli. These count rates are constant with radius within 
statistics for R = to 5 , which means that there is no 
evidence for a centrally brightened source, and hence no 
evidence for a plerionic structure centered on the pulsar. 

3. Upper limits 

The 3<j upper limits for plerionic emission confined in cir- 
cles, as well as annuli, are shown in Fig. ||. These upper 
limits were converted to absolute intensities, Sy , using the 
observed counts from Star 1 (see Fig. |l|) of Chakrabarty 
& Kaspi (1998) with V = 17 A mag. 
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Fig. 3. The 3<r upper limit to the optical flux, Sy (mJy), 
vs radius, R (arcminutes) from the radio pulsar position. 
The 'filled triangles' and the 'stars' represent annuli from 
1 and 2 arcminutes, respectively. The 'open squares' rep- 
resent the enclosed region inside each circle with radius 
R. 



The 3cr Sy upper limits increase with R as expected, 
and at R = 0. 5 (compact X-ray nebula), our 3<r upper 
limit is V > 15, while for an extended nebula of about 
4 arcminutes, the upper limit is V > 12. 
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Fig. 4. Plots of magnetic field, B (/xG), against radial dis- 
tance, R (arcminutes) from the pulsar. The thick solid line 
(lower limit) represents the pulsar driven field strength, 
assuming a a = 1 unshocked wind (Kennel & Coroniti 
1984) at large radii. For comparison, the Broadbent et al. 
(1990) and Heiles (1996) Galactic magnetic field strength 
values are also given. 

Assuming that the synchrotron spectral index, a, of 
the nebula is the same in both the optical and X-ray re- 
gions, the monochromatic X-ray flux, SikeV , at 1 keV 
is given by SikeV = Sv(e x /e )~ a , where e x (= 1 keV) 
and e (= 2.3 eV) are the X-ray and optical energies, re- 
spectively. The observed optical flux can be expressed as 
S(V) = Sy exp(— r v \ where r v is the optical depth due 
to interstellar absorption, with r v ~ 1.8 for d = 1.8 kpc. 
Using the above relations and the de Jager et al. (1995) 
expression for the TeV 7-ray IC energy spectrum (their 
Equation 2), we calculate the upper limits to the mag- 
netic fields for different values of a between 1.0 and 2.5, 
inclusively. Fig. ^ shows the curves for a = 1.0, 1.5, 2.0 
and 2.5, as well as the magnetic field limits from the X- 
ray upper limits for R = 2' for similar values of a. The 
magnetohydrodynamic curve of Kennel & Coroniti (1984) 
model, with <j = 1.0 (equipartition between particles and 
fields), is superimposed on Fig. ||, and it is considered 
as the lower limit to the extended nebula at large radii 
(R > 2/0). Also shown in Fig. || are the Galactic field 
strengths at the source according to Heiles (1996) and 
Broadbent et al. (1990). 

4. Discussions 

Our search for an extended nebula synchrotron counter- 
part to the TeV nebula did not reveal any nebular struc- 
ture around PSR B1706-44. The TeV observations do hint 
at the existence of such a pier ionic structure, but our re- 
sults constrain the nebular magnetic field strength if we 



assume that the TeV 7-rays originate from the IC scatter- 
ing of energetic pulsar wind electrons on the CMBR. At 
R = 2', we find that the optical and X-ray limits are com- 
plementary: the spectral index is constrained to a = 1.5 
to 2.3, whereas B < 3 jaG. This limit is barely above the 
pulsar wind solution, but below the average Galactic field 
values shown in Figure 4. Furthermore, our constraints on 
a are marginally consistent with a = 1.3 ± 0.3 measured 
for the compact X-ray nebula (Finley et al. 1998). 

Our optical observations and the information gathered 
from other wave bands (mainly from X-ray and 7-ray en- 
ergies) show that better quality observations (dark moon) 
should improve these limits, or, result in a detection. Po- 
larimetry of the "extended nebula" of PSR B 1706-44 may 
also help to clarify the status of synchrotron emission. 
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